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ABSTRACT 

The solar wind is connected to the Sun's atmosphere by flux tubes that are rooted in an ever-changing pattern 
of positive and negative magnetic polarities on the surface. Observations indicate that the magnetic field is 
filamentary and intermittent across a wide range of spatial scales. However, we do not know to what extent the 
complex flux tube topology seen near the Sun survives as the wind expands into interplanetary space. In order to 
study the possible long-distance connections between the corona and the heliosphere, we developed new models 
of turbulence-driven solar wind acceleration along empirically constrained field lines. We used a potential-field 
model of the Quiet Sun to trace field lines into the ecliptic plane with unprecedented spatial resolution at their 
footpoints. For each flux tube, a one-dimensional model was created with an existing wave/turbulence code that 
solves equations of mass, momentum, and energy conservation from the photosphere to 4 AU. To take account 
of stream-stream interactions between flux tubes, we used those models as inner boundary conditions for a 
time-steady MHD description of radial and longitudinal structure in the ecliptic. Corotating stream interactions 
smear out much of the smallest-scale variability, making it difficult to see how individual flux tubes on granular 
or supergranular scales can survive out to 1 AU. However, our models help clarify the level of "background" 
variability with which waves and turbulent eddies should be expected to interact. Also, the modeled fluctuations 
in magnetic field magnitude were seen to match measured power spectra quite well. 

Keywords: interplanetary medium - magnetohydrodynamics (MHD) - solar wind - Sun: corona - Sun: mag- 
netic topology - turbulence 



1. INTRODUCTION 

The Sun is known to vary on scales that span at least 
twenty orders of m agnitude in time— from milliseconds 
(Basti an et al.l [l998) to gigayears (Giidel 2007). The corre- 
sponding spatial variability may not have such a huge dy- 
namic range, but even the best observations of the solar disk — 
now reaching below 10~ 4 times the solar radius (Rq) — have 
not yet fully resolved the existing magnetic structures. There 
is also complex variability in the particles and electromag- 
netic fields measured far from the Sun in interpl anetary space 
dBruno & Carbonel 12005: Petrosya n et al.ll2010h . It is a ma- 
jor goal of solar and heliospheric physics to better understand 
how the changing conditions on the solar surface give rise to 
variations at greater distances. 

It has been suspected for many years that some trace of 
the filamentary flux-tube structure of the solar corona — e.g., 
polar plumes, jets, streamer cusps — may su r vive into in- 
terpla n etary space dMcCracken & Nessl 1 19661 iThieme et al.l 
1990L iNeu gebauer et al.1 11995b iReisenfeld et al l 119991 
Yamauchi et al. 2002; DeForest et al. 2005; iGiacalone et al.1 
2006; Borovsky 200]| 120101) . However, since heliospheric 
turbulence appears to be in a reasonably developed state, it 
is unclear how pristine solar flux tubes can avoid being com- 
pletely dispersed or broken up by the chaotic turbulent eddies 
ilMullanl [19911 iMatthaeus etaT1 ll995). Some simulations of 
magnetohydrodynamic (MHD) turbulence have shown the 
ability to prese rve input signals at single frequencies (e.g., 
Ghosh al. 2009), so it is worthwhile to explore to what extent 
the "signals" of solar flux tubes may survive as well. 

The generation of small-scale structures in the Sun's mag- 
netic field depends to some degree on the physical pro- 
cesses that heat the corona and accelerate the solar wind. 
Although these problems still have no universally accepted 
solutions, there are two broad types of explanation be- 



ing debated. First, there are models that rely on waves 
and turbulent motions to propagate up the open-field flux 
tubes and dissipate to generate the required thermal energy 
(IMatthaeus et alj|1999t ICranmer et all 120071: IChandran et all 
1201 U iVerdini et al.M2012ai) ~ Second, there are models that 
assume the wind's mass and energy is injec ted from closed- 
field region s via magnetic reconnection dFisk et al.l 119991; 
iMoore et akl 120111; lAntiochos et al.l 1201 11) . Regardless of 
whether the dominant coronal fluctuations are wave-like or 
reconnection-driven, they appear to be generated at small spa- 
tial scales in the lower atmosphere and are magnified and 
"stretched" as they evolve outward into the heliosphere. Their 
impact on the solar wind's energy budget depends on the 
multi-scale topological structure of the Sun's magnetic field. 

In this paper, we aim to improve our understanding of 
the origin of the turbulent solar wind by studying the com- 
plex structure of open coronal flux tubes. We will use high- 
resolution multidimensional models of solar wind accelera- 
tion to begin addressing several of the following questions. 
For example, how much spatial resolution is really required 
when modeling the field lines that connect the solar surface 
and the solar wind? Are there specific in situ MHD fluctu- 
ations that can be attributed to the survival of coronal flux 
tubes? Does the production of corotating interact ion regions 
(CIRs) act as a net source for turbulence dZank et al.l 1 19961) 
or is it mainly a way of smearing out l ongitudinal gradients 
and thus reducing the var iability (e.g., Smith & Wolf el [19761: 
iGoslin g 1996; Rilev 2007)? Lastly, does the well-known em- 
pirical corr elation between wind sp eed and magnetic flux tube 
expansion (Wang & Sheelev 1990) remain valid for topologi- 
cally complex bundles of field lines rooted in a mixed-polarity 
"magnetic carpet?" 

In Section[2]of this paper, we present a high-resolution po- 
tential extrapolation of the magnetic field from a quiescent 
source region of slow solar wind. We also summarize the ra- 
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dial and longitudinal properties of the open flux tubes in this 
model and discuss the solar wind speeds that would be ex- 
pected at 1 AU on the basis of existing empirical correlations. 
Section [3] describes a set of non-potential enhancements that 
we make to the field lines in order to more accurately sim- 
ulate solar wind acceleration regions. Section |4] summarizes 
the results of running these flux tubes through the turbulence- 
dr iven model of coronal heating and solar wind acceleration 
of ICranmer et al.l (12007b . In Section we describe how we 
use the plasma parameters from the one-dimensional solar 
wind models as inputs to a two-dimensional description of 
CIR formation in the ecliptic plane. We also discuss the sta- 
tistical properties of the modeled MHD fluctuations in Section 
|6l Lastly, Section [7] concludes this paper with a summary of 
the major results, a discussion of some of the wider implica- 
tions of this work, and suggestions for future improvements. 

2. HIGH-RESOLUTION QUIET-SUN MAGNETIC FIELD 

Our goal is to evaluate the importance of resolving small- 
scale magnetic features when tracing flux tubes connected to 
the low-latitude slow solar wind. Thus, we chose a time pe- 
riod during which the footpoints linked to the ecliptic plane 
were rooted in regions of Quiet Sun (QS) away from both 
unipolar coronal holes and strong-field active regions. Al- 
though QS regions are typically associated with closed loop- 
like magnetic fields, it is likely that some fraction of quies- 
cent slow wind is asso ciated with them as well (e.g. jHe et alj 
120071; iTian et alj|201 lb . Magnetically balanced QS regions, 
with only a small fraction of their magnetic flux open to the 
heliosphere, also exhibit large superradial e xpansion factors 
that a re associated with slow wind speeds (Wan g~& Sheelevl 
fl990h . 

We used magnetic field measurements made by the Vec- 
tor SpectroMagnetograph (VSM) instrument of the Synoptic 
Opt ical Long-term Investigations of the S un (SOLIS) facil- 
ity (EdiEiLiD HH Einniiitin HHI). SOLIS magne- 
tograms may not have the highest spatial resolution in com- 
parison to other available data sets, but its sensitivity to weak 
fields made it the optimal choice for mapping accurate foot- 
points of open QS flux tubes. The primary measurement was 
a full-disk longitudinal magnetogram taken in the Fe I 6301.5 
A line at time t Q = 16:30 UT on 2003 September 4. For global 
context, the central portions of this magnetogram were em- 
bedded in a lower-resolution synoptic magnetogram for Car- 
rington Rotation (CR) 2007. At time to, the central meridian 
longitude during this CR was (j>o = 286.6°. 

The magnetic field was extrapolated from the photosphere 
to the corona using the standard Potential Field Source Sur- 
face (PFSS) method, which assumes the corona is current- 
free out t o a spherical "source surface" above which the field 
is radial (ISchatten et al.lll969t lAltschuler & Newkirklfl969b . 
The source surface radius was chosen as r = 2.46/?©. De- 
tails about the numerical method used to construct the PFSS 
model are provided in Appendix [A] The PFSS method has 
been shown to create a r elatively good mapp i ng between the 
Sun and the heliosphere (Arge & Pizzo 2000; Luhman n et al.l 
12001 IWang & Sheelevl 120061) . although full MHD simula- 
tions of course take better account of gas pressure effects and 
stream-stream interactions (see Section |5). 

We traced a set of "open" field lines down from the source 
surface to the photospheric lower boundary. The locus of 
starting points at r = 2A6Rq was coincident with the eclip- 
tic plane. The longitudinal separation in the grid of starting 
points was chosen to be equivalent to 1 minute of solar rota- 




Figure 1. (a) View from above the ecliptic plane of PFSS field lines traced 
from the source surface down to the photosphere. Central regions of the high- 
resolution SOLIS magnetogram are embedded in a lower-resolution synop- 
tic magnetogram for CR 2007. (b) Saturated grayscale version of the high- 
resolution magnetogram with lines of heliographic longitude and latitude 
shown at 10° separation (yellow dotted curves). Also shown are footpoint 
locations of the 8727 high-resolution field lines (red points) and of field lines 
traced from the same starting locus, but using low-resolution (<? max = 9) WSO 
synoptic data (blue points). 

tion time. The total longitudinal extent of this region was 80°, 
which placed the footpoints squarely inside the disk-center re- 
gions of the high-resolution magnetogram obtained on 2003 
September 4. Using a Carrington rotation period of 27.2753 
days to define the angular rotation rate (il = 2.6662 x 10~ 6 
rad s -1 ), we produced a total of 8727 field lines separated in 
azimuthal angle by Ac/) = 0.00917°. 

Figure [TJ shows two illustrative views of the field lines and 
footpoints. A subset of the field lines is shown from a view- 
point above the ecliptic plane in Figure 01 a X along with the 
high-resolution SOLIS magnetogram embedded in the low- 
resolution synoptic magnetogram. Figure [Ub) shows the full 
set of 8727 footpoints on a saturated version of the high- 
resolution magnetogram (i.e., three gray levels only) from the 
viewpoint of an observer on Earth. We also plot an alternate 
locus of footpoints that was computed from a low-resolution 
PFSS model made with synoptic mag netogram data from the 
Wilco x Solar Observatory (WSO; see Hoeksema & Scherrerl 
119861) . The published PFSS coefficients for CR 2007 were 
constructed with maximum order £ = 9 in the spherical har- 
monic expansion. FigureQJb) shows that at this time the heli- 
ographic latitude Bo was near its annual maximum of +7.2°. 

The magnetic polarities in the QS footpoint area of in- 
terest were largely mixed and balanced, but there was a 
slight preponderance of negative polarity in this region of 
the photosphere. An approximate measurement of the flux 
imbalance fraction £ (i.e., the ratio of net flux density to 
absolute unsigned flux density) yielded £ « 0.15, which is 
well within the range of values expected for QS regions 
( Wieg elmann & Solankil 12004 lHagenaar et a l. 2008). The 
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slight imbalance toward negative polarities was amplified 
at larger distances, such that at the PFSS radius the entire 
80° wide sector shown in Figure 02a) was uniformly nega- 
tive in its polarity. Thus, each of the 8727 open field lines 
exhibited a negative polarity at their footpoints as well. 

If the full set of field lines was mapped down radially from 
their initial locations to the solar surface, the horizontal sep- 
aration of each neighboring pair of footpoints would be ap- 
proximately 110 km. It is clear from Figure Q~]that the actual 
distribution of all 8727 footpoint separations is quite broad. 
Clumps of closely packed footpoints are separated by wider 
jumps, the latter occurring in regions where the large-scale 
polarity changes. The distribution of neighboring footpoint 
separations spans several orders of magnitude, with mini- 
mum, mean, and maximum values of 5.2, 281, and 249,000 
km, respectively. The distribution has a median value of 40 
km, and it is skewed such that 75% of the neighbor separa- 
tions are less than 100 km. 

In contrast to the statistics of neighboring footpoint sep- 
arations given above, we note that the SOLIS pixel size is 
roughly 820 km on the solar surface. Thus, our reconstructed 
magnetic field lines often tend to oversample the information 
available in the photospheric data. However, there is evi- 
dence that the overlying coronal magnetic field may exhibit 
topological features — such as nulls, separatrices, and quasi- 
separatrix layers (QSLs) — that exist on s patial scales smaller 
than those of the driving s ources (e.g.. lCloseeTaT1 l2005b: 
lAntiochos et alJl2007L 1201 ll) . It is not clear how much extra 
resolution is needed in order to resolve these features prop- 
erly and to produce models that contain the full diversity of 
flux-tube expansion properties in the open coronal fie ld, so 
we chose to over-resolve (see also Garraffo et al. 2012j). 

Figure [2] shows the radial dependence of magnetic field 
strength along a subset of the high-resolution PFSS field lines. 
Above the source surface radius of 2.46 R®, the radial field 
strength is extrapolated as B r oc r~ 2 . The models described in 
this section contain no explicit azimuthal field component 
above the source surface. At 1 AU, the modeled field strengths 
range between about 10~ 5 and 2.5 x 10~ 5 G, which is about a 
factor of two smaller than i s generally measured in the ecliptic 
(Maria ni & Neub auer 1990). However, the two-dimensional 
models described in Section contain both a self-consistent 
Parker spiral field (i.e., B^ ^0) and intermittent enhancements 
of B, due to stream interactions. 

For the full set of 8727 modeled field lines, the distribution 
of field strengths at the photospheric lower boundary has a 
minimum and maximum of 1 . 1 6 and 1 3 1 G, respectively. The 
median photospheric field strength is 24. 1 G, the mean is 3 1 .3 
G, and the standard deviation about the mean is 26.3 G. Fig- 
ure |2] shows 36 field lines from the PFSS model: 34 chosen 
at random, plus the ones with minimum and maximum field 
strengths at the photosphere. In Figure [2] we also plotted the 
radial dependence of field strength s for t he flu x tubes used in 
the solar wind models of Cranmer et al. (20071). 

We use the modeled radial distributions of field strength to 
extract a convenient scalar measure of each flux tube's cross- 
sectional expansion in the corona. In order to quantify the 
degree of superradial expansion, we define the Wang-Sheeley- 
Arge (WS A) expansion factor / ss as 



where the subscript "base" refers to the photospheric lower 
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Figure 2. Radial dependence of magnetic field strengths along a subset of 
field lines traced in the SOLIS PFSS model (black s olid curves). We comp are 
these with the magnetic field distributions used by Cranmer et al. 1 2007) for 
ZEPHYR solar wind models. Solid red curves show a range of latitu des in the 
quiet solar minimum configuration of Banaszkiewicz et al. ( 1998). Dashed 
red curves show open flux tubes assumed to be rooted in strong active regions. 
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Figure 3. Empirical relationships between the flux-tube expansion fac- 
tor / ss and the terminal sol ar wind speed Moo. Publis hed parameteriza- 
tions of Arge & Pizzo 12000) (red dot-dashed curve) and Riley & Luhmann 
120121) ( blue dashed curv e ) are compared with the original data-ranges 
given by Wane & Sheelev (1990) (gray region), solar wind model results 
of Cranmer et al. (2007) (green dotted curve), and Equation l[2} (black solid 
curve). 

boundary of the PFSS extrapolation and "ss" refers to the 
source surface radius of 2.46 Rq. Much like the footpoint lo- 
cations and photospheric field strengths, we find that / ss varies 
intermittently with longitude. 

There are several independent calibrations of the well- 
known empi rical relationship between f 5S and the sola r 
wind speed dWang & Sheelevl 119901: lArge & Pizzol 120001) . 
These parameterizations depend on the spatial resolution 
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and noise level of the photospheric magnetogram data, as 
well as on the method used to extrapola t e the field (see 
also [Wang etaTI 119971: lOwens et all [20051; iMcGregor etaTI 
1201 11 iRiley & Luhmannll2012l) . Figure [3] compares several 
of these parameterization s with one another. Note that the 
IRiley & Luhmannl (120121) curve largely repro duces the range 
of values r eported by the [Wang et all (11997b revaluation of 
the original IWang & Sheelevl (119901) analysis. Figure [3] also 
shows a "concordance" WS A relationship that we constructed 
to fall within the range of the existing correlations. This rela- 
tionship is given by 



o 100 - 



2300 



ln/ ss + 1.97 



km s 1 



(2) 



where is the terminal or asymptotic solar wind speed. We 
do not take into account the other commonly used variable 
parameter of the WSA model: the transverse angular distance 
9b between the field line of interest and the boundary of the 
nearest large-scale coronal hole. When using this parameter, 
it is often found that the empirical wi nd speed becomes inde - 
pendent of this angle for > 7° (e.g.. lMcGregor et"aDl20 1 ll) . 
The patch of QS in which we model the footpoints of slow so- 
lar wind appears to be far enough away from any large coronal 
holes that these field lines are probably insensitive to Ob- 

Figure|4]shows the longitudinal dependence of the flux-tube 
expansion factor and empirical wind speed for all 8727 
modeled field lines. We also show corresponding parameters 
of the PFSS model computed from the low-resolution WSO 
magnetograms. Figure |4ja) plots / ss as a function of a lon- 
gitudinal angle (0 - <po) normalized by the central-meridian 
Carrington longitude from the high-resolution SOLIS mea- 
surement of 2003 September 4. Figure @Ib) shows Uoq, com- 
puted from Equation d2j, as a function of time t estimated for 
solar wind flux tubes to rotate past the Earth. For this plot, the 
conversion between longitude and time is given by 



-ra 



?o+- 



Ar 



(3) 



where fn is the time that longitude 0o was on the central merid- 
ian, and we use Ar = 214^ Q and (u) = 450 km s" 1 to estimate 
the mean travel-time from the Sun to 1 AU. Time is plotted 
in Figure 2] in units of the day of year (DOY) in 2003, and 
the modeled 80° of longitude converts into approximately 6 
days of rotation time. In Figure Hfc) we show a subset of 6 
hours to illustrate that extremely sharp structures in / ss exist 
in the PFSS model that may generate similarly sharp struc- 
tures in the solar wind (at least close to the Sun). In Section 
[5] we model the transit-time evolution of the solar wind more 
accurately. 

We emphasize that our high-resolution PFSS extrapola- 
tion of the Sun's magnetic field represents merely a single 
"snaps hot" in time. In reality , the small-scale magnetic car- 
pet (Title and Schriiver 1998) evolves over a wide range of 
timescales and causes the open-field footpoint locations to 
move around. Specifically, observations have indicated that 
the "recycling times" for magnetic flux in the phot osphere and 
corona probably are between 1 and 10 hours (e.g. JClose et al.l 
I2005at lHagenaar et all 12008b . Monte Carlo models of the 
magnetic carpet's evolution also found a similar range for the 
mean times cale of the opening up of closed flux tubes via re- 
connection (Cranmer & van Ballegooijen 2010). These times 
are comparable with how long it takes the radial projection of 
a supergranular cell to rotate past a distant observer (roughly 
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Figure 4. (a) Longitudinal dependence of / ss computed from the high- 
resolution SOLIS model (black solid curve) and the low-resolution WSO 
model (red solid curve), (b) Empirically estimated wind speed «oo shown 
as a function of time in day-of-year (DOY) units, with same curve types as 
in (a), (c) Six-hour expanded subset of the data shown in (b). The central 
meridian location of the SOLIS high-resolution magnetogram is shown in (a) 
and (b) with a blue dotted line. 

3^1 hr). Thus, we do not mean to present the high-resolution 
PFSS reconstruction as a true dynamical model, but only as 
a representative state of the QS field that the solar wind will 
"see" as it accelerates up through the open flux tubes. 

3. NON-POTENTIAL RADIAL ENHANCEMENTS TO THE FIELD 

The previous section described the straightforward PFSS 
extrapolation of measured photospheric fields into the corona. 
However, that process does not take into account the full range 
of magnetic field variations that we believe exist along these 
flux tubes. In this section we describe several adjustments 
that are made to the radial dependence of the magnetic field 
strength B{r) for the modeled set of time-steady field lines. 

First, we recognize that the photospheric footpoints of the 
large-scale coronal magnetic field appear to be broken up into 
thin flux tubes (i.e., observed widths of order 50-200 km) 
that collect in the dark lanes between the ^1000 km diam- 
eter granulation cells. These flux tubes have field strengths of 
1-2 kG and are often called "G-band bright points" (GBPs) 
because they show up as bright features in the 4290-4320 A 
molecular bands (e.g., Bergeret al. 19951 IS teiner et al.1l2001h . 
Horizontal motions of these features have been used to put 
empirical constraints on the photospheric flux of Alfven wave 
energy that propagates up into the corona dNisenson et al.l 
l2003TlCranmer & van Ballegooiienll2005l: IChitta et al.ll2012lh 
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Figure 5. Schematic illustration of flux-tube expansion in the presence 
of non-potential fields. In both panels, open flux tube boundaries (blue 
solid curves) surround empirically determined field-free regions (gold ar- 
eas), (a) For Bi ^0 and Bi = 0, the only enhancements to the open field 
occur below the mid-chromosp heric "merging height" (dashed curve) of 
Cranmer & van Ballegooijen (2005). (b) Postulated field enhancements in 
both the photosphere and the low corona are shown (i.e., Bj /0 and Bi ^0). 

GBPs do not show up individually in the SOLIS magne- 
tograms that we used to reconstruct the coronal field, so we 
modify B(r) as described below to account for their supposed 
presence. 

Figure |5fa) illustrates how the assumption of photospheric 
field fragmentation results in narrower flux tubes and higher 
field strengths than would be otherwise obtained from the 
magnetogram data. Somewhat crudely, we assume the lower 
solar atmosphere is divided into strong-field flux tubes and 
"field free" regions that are in total pressure equilibrium with 
one another. Enhancing the field strength inside the flux tube 
is equivalent to assuming a larger field-free volume between 
the tubes. Thus, along each flux tube we modify the origi- 
nal PFSS field Bo(r) by adding two additional components in 
quadrature; i.e., 

B 2 = B\ + B\ + B\. (4) 
The photospheric GBP enhancement B\ is given by 

B 1 (z) = B Q exp(-^j , (5) 

where z = r—R® is the height above the (optical depth unity) 
photosphere, and we adopt Bq = 1470 G as a universal GBP 
field strength. The upper photospheric scale height Hi = 
120 km corresponds to a hydrostatic temperature of approxi- 
mately 4000 K. These constants result in a field strength that 
ma tches the photospheric and lower ch romospheric parts of 
the Cranmer & van Ballegooijen (2005) flux tube model. 

There is evidence that additional modifications to a poten- 
tial field (i.e., the B-> term above) are needed in the corona. For 
example, Cran mer & van Ballegooijen (2005) took account of 
the large-scale bundling of fields in the supergranular network 
of a coronal hole by constructing a two-dimensional magne- 
tostatic mo del of the canopy-like expansion of open flux tubes 
(see a lso iGabriell [T976t iGiovanellil [l980t lHackenberg eTaTI 
20001 lAlouaz & Rasj 1200^ 7 The presence of weak-field re- 



gions in the internetwork cell centers (with presumably higher 
gas pressure than the strong-field network lanes) modifies the 
potent ial field expansion at chromospheri c heights. Simi- 
larly, ISchrijver & van Ballegooijenl (120051) studied the modi- 
fications to a potential field due to the presence of finite gas 
pressure gradients throughout the low corona. They modeled 
the spatial dependence of the plasma /3 parameter (i.e., the ra- 
tio of gas pressure to magnetic pressure) in a QS region and 
found that ^fsl over much of the volume. This implies the 
gas influences the field topology in ways unanticipated by the 
potential field model, which implicitly assumes j3 <C 1. 

It is well known that active regions frequently contain large- 
scale currents that give rise to twisted , sigmoidal field lines 
(g-g^ lGarv etail H987t iCanfield et al.1 [T999t iSchriiveret all 
2005). These structures are believed to be the result of a 
combination of surface shear motions and the emergence of 
new flux from below the solar surface. However, because 
these effects also occur elsewhere on the Sun, it is likely that 
other regions (including the QS footpoints of open flux tubes) 
exhibit currents a nd non-potential fields on a wide range of 
spatial scales (see[Ab bett 2007 ;IZhao et al.l2009HYeates et al.l 
I2010t iReardon et alj|201 U iMever et al.ll2012l) . These effects 
give rise to increased "fibril" type complexity to the field. 
Their presence is likely to increase the field strength in the 
low corona and thus shrink the volume of any given open flux 
tube that traverses the non-potential region (see Figure 0b)). 
Other physical effects that may contribute to modifying the 
radial dependen ce of B include per vasive chromospheric up- 
flows (Mcintosh & De Pontieu 2009), rotational supergranule 
motions (Zhang & Li ul201 11) . and loop footp oint asymmetries 
that g ive rise to rapid unresolved motions (IWiegelmann et all 
I20T1 . 

In the absence of a clear-cut method of modeling non- 
potential field enhancements in the low corona, we adopt a 
similar hydrostatic radial dependence as was used in Equa- 
tion (|5]l. The added magnetic field component is given by 

B 2 (z) = B Q (0)exp(-^j , (6) 

where an approximate scale height H2 = 2k^T2/{m p g) is de- 
fined by an arbitrary effective temperature T2, the Boltz- 
mann constant k^, the proton mass m p , and the Sun's sur- 
face gravitational acceleration g[3 Because the primary goal 
of this non-potential enhancement is to "stretch out" the pre- 
existing field, it is normalized using the PFSS lower bound- 
ary condition Z?o(0) of that particular field line. The ef- 
fective temperature T2 is a purely empirical parameter that 
characterizes the radial extent of the B2 field modification. 
The only way that 7i would be related to an actual coro- 
nal temperature would be if gas pressure effects were the 
primary cause of the non-potentia l enhancement (see, e.g., 
ISchrrjver & van Ba llegooijen 20051). 

Figure |6j a) illustrates the magnetic field enhancements de- 
scribed above using a mean PFSS model. The unmodified 
field strength Bo, shown by the dashed curve, was produced 
by forming the average of In Bq at each height and taking the 
exponential of the result. The photospheric enhancement B\ 
dominates the modified field strength at heights z < O.OO2/? . 
The other curves show the result of varying the effective scale- 
height temperature T2 in Equation ©. Values of T2 between 

1 We simplified the definition of H2 by assuming a fully ionized hydrogen 
plasma; i.e., we did not include the impact of helium and other heavy ions on 
the mean molecular weight, and we did not include partial ionization effects. 
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Figure 6. Radial dependence of (a) mean magnetic field strengths and (b) 
ratios of rms field strength fluctuations to the mean field strengths. The un- 
modified Bo model (black dashed curve) is shown with a model that includes 
only the B\ enhancement (black solid curve) and others that also contain the 
Bo coronal enhancement (color solid curves, see labels for T% values). Also 
shown in (b) is the ratio of rms to mean field strength in the non-WKB Alfven 
wave model of Cranmer & van Ballesooiien (2005) (gray dotted curve). 

about 0.05 and 0.5 MK represent mild enhancements to the 
field strength in the low corona that are consistent with the 
flux-tube constrictions illustrated in Figure HJb). The red 
curve corresponding to the largest field strength in Figure |6] 
which was computed for T2 = 1 .5 MK, is probably beyond the 
realm of physical realism. Generally, we do not expect the 
non-potential effects described above to extend as far up as 
z > 0.5^0 in a QS region. We also do not believe that using 
a value of T2 that is of the same order of magnitude as the 
expected peak coronal temperature is especially realistic, ei- 
ther. We show the Tz = 1.5 MK curve mainly for the sake of 
completeness. 

In Figure |6jb) we show a measure of the statistical vari- 
ability of the field strength in the whole set of 8727 SOLIS 
flux tubes. The plotted quantity is the ratio of the standard 
deviation in B to its mean value at a given height. As above, 
the dashed line shows the unmodified field strength Bq, which 
reaches a maximum ratio of ^0.84 in the photosphere and de- 
clines to about 0.27 at the source surface. When the B\ and Bo 
enhancements are added there is necessarily a smaller degree 
of variability in B, especially for the strong photospheric com- 
ponent at low heights. In the corona, the ratio gets smaller for 
two reasons of nearly comparable magnitude to one another: 
(1) the mean field strength increases, and (2) the standard 
deviation decreases. We also compare these ratios of rela- 
tive variabil ity to a representative model of Alfv en wave vari- 
ability from Cranme r & van Ba llegooijen (2005). The plotted 
quantity is the ratio of the root mean squared (rms) Alfvenic 



fluctuation amplitude SB to the background field strength Bo. 
Of course, this is somewhat an "apples to oranges" compar- 
ison because the Alfvenic variability represents propagating 
transverse waves and the other curves describe a collection of 
static spatial fluctuations. Nonetheless the overall similarity 
of the two ratios indicates that the two types of variations may 
be of comparable importance in the low corona. 

We chose not to include the non-potential field enhance- 
ments B\ and B2 in our calculations of the superradial expan- 
sion factor / ss for each field line. These enhancements are 
included mainly for the benefit of the ZEPHYR simulations, 
which make use of the Alfven speed profile in the lower atmo- 
sphere to translate photospheric velocity amplitudes into wave 
energy fluxes. If we had inserted the modified field strengths 
in Equation (Q~|i to rescale / ss , the basal field strengths would 
all be nearly identical to one another due to the B\ term. In 
that case, much of the relevant information about field-line 
spreading in the upper chromosphere and low corona (which 
is contained in the measured variability of Bo(0)) would have 
been erased from / ss . 

4. TURBULENCE-DRIVEN SOLAR WIND MODELS 

We used the computed magnetic fields B{r) as inputs to a se- 
ries of one-dimensional physical models of turbulent coronal 
heating and solar wind acceleration. The steady-state models 
presented below are numerical solutions to one-fluid conser- 
vation equations for mass, momentum, bulk inter nal energy, 
and Alfven wave energy. ICranmer et al.l ([2007 ) described 
these equations and outlined the computational methods used 
to solve them using a computer code called ZEPHYR. The 
models presented below were calculated with a slightly mod- 
ified version of the original ZEPHYR code. In addition to 
several algorithmic improvements that were needed to allow 
the code to read in large numbers of flux tubes at a time, we 
made three modifications that affected the results: 

1 . We chang e d the value o f the coefficient that multiplies 
iHolrwegT s dl974fl9 76) prescription for free-streaming 
conductive he a t flux in the collisionless heliosphere. 
ICranmer et al.1 (120071) used a c = 4, as originally sug- 
gested by Hollweg, but we reduced it to a c = 1 based 
on a recent analysis of ele ctron heat flux measu rements 
in the fast solar wind (see ICranmer et al]|2009l) . 

2. We reduced the photospheric boundary condition on the 
energy flux of longitudinal acoustic waves from 10 8 to 
10 6 erg s -1 cm" 2 . The lower value gave a more real- 
istic height for the transition region (TR) between the 
chromosph ere and corona than did the higher value. 
Figure 8 of Cr anmer et al.l (120071) showed that a larger 
acoustic wave pressure in the chromosphere gives rise 
to a larger density scale height and thus causes the "crit- 
ical" density for runaway radiative instability to occur 
at a larger height. Our adopted value of the acous- 
tic wa ve flux, in combination with the Cranmer et al. 
(2007) choice for the photospheric Alfven wave ampli- 
tude (vj_ = 0.255 km s" 1 ), was held fixed for all of the 
models discussed below. 

3. We used the modified version of the numerical relax- 
ati on method for th e internal energy equation described 
by ICranmer] (12008). We also reduced the initial value 
of the minimum undercorrection exponent eo from 0. 17 
to 0.10. These changes gave rise to more robust conver- 
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gence of the coronal and heliospheric temperature T(r) 
to its steady-state solution. 

It should also be emphasized that ZEPHYR code makes use 
of the full radial dependence of B(r) and does not depend on 
the spatial resolution of observations that can affect the nor- 
malization of the expansion factor / ss . 

Prior to computing models for a large number of flux 
tubes from the SOLIS reconstruction, we performed a lim- 
ited parameter study to explore the effects of varying the non- 
potential field enhancement Bi- We began with the mean-field 
models shown in Figure |6ja) and produced a finer grid that 
varied the T2 parameter from to 1.5 MK in increments of 
0.025 MK. These models all used identical lower boundary 
conditions in the photosphere, and only differed in their tabu- 
lated B(r) field strengths. Of those 61 models, 52 of them con- 
verged successfully to a steady-state solution that satisfied in- 
ternal energy conservation to within 5% accuracy Q The other 
9 models corresponded to values of Ti > 0.9 MK, which we 
suspect is outside the realm of physical realism for the non- 
potential field enhancements. 

In Figure [7] we show several summary parameters of the 
successful subset of mean-field ZEPHYR models that varied 
the T2 parameter. As the non-potential field strength in the 
low corona is increased, the wind speed u at 1 AU decreases 
and the proton number density n p at 1 AU increases. Fig- 
ure [7Ja) shows that most of this variation occurs as T2 is in- 
creased from to 0.5 MK. Subsequent increases from 0.5 to 
1 .5 MK do not appear to produce substantial changes in the 
modeled solar wind. Figure |7fb) plots the maximum coro- 
nal temperature versus the Alfven wave velocity amplitude 
v_i_ at 1 AU. One can think of this value of as a "resid- 
ual" amplitude since it is the end result of wave dissipation 
that occurred in the corona and inner heliosphere. It makes 
sense that larger coronal temperatures correspond to smaller 
values of v_i_, since more coronal heating is consistent with 
more damping of MHD turbulence. 

We can explain the variations shown in Figure|7]by noticing 
that the turbulent heating rate Q is often close to being linearly 
proportional to the background field strength B in the upper 
chromosphere and low corona. ICranmerl (t2009l) showed that 
this proportionality is exact in cases where both the thin flux- 
tube relation (B oc p l l 2 ) and sub-Alfvenic wave action conser- 
vation (vj_ oc p -1 / 4 ) are valid. Thus, when the field strength 
Z?2 is increased, the amount of heat deposited at the coronal 
base increases as well. In direct response, the corona's base 
pressure increases, as does the solar wind's mass loss ra te (see 
lHammeri[T98l lWimbrod[l988l iCranmer & Saa^l20TTl) . This 
explains why the peak coronal temperature and the proton 
number density at 1 AU both increase when T2 is increased. 
The decrease in v± at 1 AU was explained above as a result 
of the increased wave damping that goes along with stronger 
coronal heating. 

The variation in the wind speed u in Figure [7] can be un- 
derstood as a result of th e presence of p onderomotiv e wave - 
pressure acceleration (Jacquei 119771: iHolzer et al.l 11983b . 
When the relative strength of coronal heating is low (i.e., for 
the smallest values of T2), the wind is driven primarily by 
wave pressure and not gas pressure. In that regime, there is 
a higher outflow speed when there is a stronger (less damped) 

2 In other words, these models exhibited final values of the convergence 
p arameter (8E) less tha n 0.05. This parameter was defined in Equation (63) 
of Cranmer et al. (2007) and its iterative convergence was illustrated in Fig- 
ure 4 of that paper. 
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Figure 7. Results of mean-field ZEPHYR solar wind models for a range 
of T2 parameters, (a) Wind speed at 1 AU in units of 100 km s (solid 
symbols) and proton number density at 1 AU (open symbols). Symbol colors 
are mapped directly to the T2 values shown on the horizontal axis, (b) Alfven 
wave velocity amplitude at 1 AU plotted vs. the peak coronal temperature for 
the mean-field models (symbol s with same T2 colors as in panel (a)) and for 
the models presented by Cranmer et al. (2007) (gray region). 

population of Alfven waves both in the corona and at 1 AU. 
Figure |7Jb) shows that when T2 is smaller than about 0.25 
MK, the corona is heated to peak temperatures below 1 MK 
and the res idual wave amplitud e at 1 AU is higher than was 
seen in the lCranmer et al.l (120071) ZEPHYR models. Thus, we 
may be able to rule out values of T2 below ^0.25 MK because 
of their unrealistic solar wind properties. If we combine this 
with the discussion above that appeared to also rule out values 
of T2 larger than 1—1.5 MK, this leaves only a limited range 
of T2 parameters (i.e., 0.25 < 7s < 0.75 MK) that may be rel- 
evant for modeling our reconstructed QS flux tubes. 

With the above results in mind, we proceeded to create so- 
lar wind models for the SOLIS flux tubes discussed in Section 
12 Because the full set of 8727 flux tubes often oversamples 
the existing complexity of the reconstructed magnetic field, 
we selected a subset of 289 field lines to input to ZEPHYR. 
These field lines were chosen by sampling the data more fre- 
quently in regions of rapid longitudinal change than in regions 
where the field strength was nearly constant. Every local max- 
imum and minimum in the longitudinal variation of / ss (see 
Figure Ufa)) was represented by at least one of the 289 re- 
sampled points. We ran ZEPHYR for three distinct cases of 
non-potential field enhancement: (1) T2 = (i.e., using the 
photospheric B\ enhancement only), (2) T2 = 0.2 MK, and (3) 
T2 = 0.5 MK. We decided against using even larger values 
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Figure 8. Comparison of asymptotic solar wind speeds estimated from 
the WSA anticorrelation (lower black curve) with those computed from 
ZEPHYR models having T2 = (upper blue curve) and T2 = 0.5 MK (middle 
red curve). Longitudes are plotted for the footpoints of the PFSS flux tubes 
as in Figure|4ja). 

of T2 because Figure [7] showed that the resulting solar wind 
would probably not have been significantly different from the 
T 2 = 0.5 MK case. 

Figure [8] summarizes the results of two of the ZEPHYR 
models (T2 = 0, 0.5 MK) by comparing the wind speeds at 
1 AU with the predicted WSA wind speed from Equation (|2). 
The intermediate case of T2 = 0.2 MK exhibits plasma prop- 
erties that lie between the two plotted models. Note that an 
overall WSA-type anticorrelation between u and / ss appears 
to be upheld for the numerical models. As was discussed 
above, the T2 = models do not contain enough turbulent 
heating to give rise to coronal temperature maxima above 1 
MK. Thus, the resulting lack of wave dissipation gives rise 
to "too much" wave pressure acceleration and a high-speed, 
low-density wind. The case of T2 = 0.5 MK gives rise to wind 
speeds at 1 AU between 448 and 853 km s" 1 , which is a bit 
higher than the observed range of speeds but is much more 
realistic than the predictions of the T2 = model. 

In Figure |9ja) we examine in more detail how the models 
appear to follow WSA-like relationships between wind speed 
and expansion factor. Although there is substantial scatter, 
each set of ZEPHYR results does seem to be centered around 
a relationship reminiscent of Equation (j2). Models with lower 
values of T2 correspond to larger normalization offsets in the 
wind speed. For the T2 = 0.5 MK case, we show the modeled 
wind speeds measured at a heliospheric distance of r = 20 Rq 
in addition to those at 1 AU. At the former distance the wind 
speed in each flux tube is roughly 75% of its asymptotic termi- 
nal speed in interplanetary space. Those speeds agree remark- 
ably well with the curve corresponding to Equation (0. This 
is consistent with the fact that the empirical WSA relation- 
ship is often applied as a lower boundary condition (typically 
at distances of or der 20 Rq) for global simulations of the inner 
heliosphere (e.g. JLee et at]l2009tlMcGregor et al.ll201 lb . 

It is difficult to pin down a simple explanation for the 
emergence of clear WSA anticorrelations in Figure HJa). 
There have been various attempts to explain how an Alfven 
wave driven wind should exhibit this kind of relationship be- 
tween wind speed and expansion fac tor (iKovalenkol Il981h 
IWang & Sheelevl I199U |Suzuki 2006). These explanations 
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Figure 9. (a) Solar wind speeds plotted vs. / ss compared with Equation (2} 
(black solid curve). Speeds at 1 AU are shown for ZEPHYR models having 
T2 = (blue filled symbols), = 0.2 MK (green open symbols), and To = 0.5 
MK (red filled symbols). Speeds at r = 20 Rq are also shown for the Ti — 0.5 
MK model (gray filled symbols), (b) Proton number densities n p at 1 AU 
plotted vs. wind speeds for the T2 = 0.5 MK model (large black symbols) and 
for AC£/SWEPAM data (small gold symbols). 



tend to involve differences in the radial evolution of wave en- 
ergy flux along different field lines — i.e., larger wave fluxes 
at the critical point correspond to both greater acceleration 
and smaller values of / ss . This must be happening to some 
extent, but the ZEPHYR models also include non-WKB re- 
flection, self-consistent wave damping, and wind accelera- 
tion from both gas pressure and wave pressure. These ef- 
fects are linked to one another via a number of nonlinear feed- 
backs, so we do not yet know the precise chain of events that 
gives rise t o the emergent distribution of wind speeds. Future 
work (e.g., Woolsev & Cranmer 2012) will aim to determine 
whether these interactions can be understood from the stand- 
point of more basic scaling relations. 

Figure |9jb) shows the anticorrelation between wind speed 
and proton number density at 1 AU. In addition to the T2 = 0.5 
MK model results, we also show hourly-averaged data mea- 
sured by the Solar Wind Electron, Proton, and Alpha Moni- 
tor (SWEPAM) instrume nt on the Advanced Co mposition Ex- 
plorer (ACE) spacecraft dMcComas et al.l!998l) . The data dis- 
played are a relatively sparse selection of approximately 4700 
data points spread out in time between the years 1998 and 
2006. Even though the models do not extend to the lowest 
wind speeds and highest densities seen in the low-latitude he- 
liosphere, the overall sense of agreement between the models 
and the measurements is good. 
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Figure 10. Radial dependence of solar wind parameters for a selection of five 
flux tubes extracted from the ZEPHYR models computed with T2 = 0.5 MK. 
Shown are: (a) one-fluid temperatures, (b) radial wind speeds, and (c) mag- 
netic field strengths vs. height above the solar photosphere. Line styles and 
colors serve merely as consistent identifying labels for the flux tubes. 

To illustrate the plasma properties that appear in the 
ZEPHYR models, we show in Figure [10] the temperatures, 
outflow speeds, and magnetic field strengths for a selection of 
five flux tubes from the T2 = 0.5 MK model. In the low corona, 
the magnetic field strengths in these flux tubes vary over about 
an an order of magnitude (from 4 to 40 G). However, Figure 
[ToT c) shows that B(r) does not vary substantially in the photo- 
sphere and chromosphere (z < 0.002 Rq) nor in the extended 
corona and heliosphere (z > 0.5 Rq) of these models. As de- 
scribed above, stronger coronal fields give rise to more heat- 
ing below the critical point, which in turn pr oduces a hotter, 
denser, and ultimately slower wind (se e also Leer & Holzer 
[19801 lPneumanl[T980l ILeer et al.lll982l) . Figure Mb) shows 
that the hotter models undergo more intense "bursts" of so- 
lar wind acceleration in the low corona, but their large-scale 
outflows end up being slower. 



Figure [Tol a) shows that the model receiving the largest de- 
gree of coronal heating has a higher TR than the other models. 
This occurs because that model's intense wave pressure in the 
upper chromosphere produces a larger pressure scale height 
and a shallower radial decline in the density. Thus, the den- 
sity does not drop to its critical value (for runaway radiative 
instability) until it reaches a significantly larger height than in 
the other models. A situation similar to the above case occurs 
for roughly half of the full set of 289 ZEPHYR models com- 
puted with Ti = 0.5 MK. Defining the TR height zjr as the 
location at which T first exceeds 0.2 MK, we found that the 
models bifurcate cleanly into two groups: 162 of them had 
Ztr < 0.008^0 and 127 of them had z T r > O.O12fl . As ex- 
pected, the latter set of models with larger values of zjr tended 
to have larger peak temperatures (r max > 1.3 MK) and lower 
asymptotic wind speeds {u^ < 600 km s" 1 ) than the other 
models. This appears to be a kind of "bistable" situation, and 
it is possible that a truly time-dependent model would exhibit 
oscillations between these two states. We will study this effect 
in future work. 

5. INTERPLANETARY EVOLUTION OF THE SOLAR WIND 

Each of the individual ZEPHYR models described above 
was computed independently of the others. However, we 
know that magnetic flux tubes in the corona and heliosphere 
interact with their "neighbors" in a variety of ways. The dis- 
tinct identity of a given flux tube can be blurred by time- 
dependent di sruptive effects such as Kelvin-Helmholtz insta- 
bilitie s (e.g.. lEinaudi"etin[T999l;[Parhi et al.lfl999i: iVelli et all 
1201 lb or diffusive field-line shredding du e to MHD turbu- 
lence dMullarill99HlMatthaeus et al.lll995l) . Also, the large- 
scale structure of the heliosphere is often dominated by colli- 
sions between obliquely oriented streams that result in CIR- 
like compressions and rarefactions. 

Our goal in this section is to begin taking account of the 
multidimensional evolution of solar wind flux tubes as ini- 
tially computed from the high-resolution magnetic field mod- 
els of Sections [2^4] We use the longitudinally resolved set of 
ZEPHYR plasma properties as an inner boundary condition 
to a two-dimensional MHD calculation in the ecliptic plane. 
Thus, we simulate the time-steady effects of CIR formation, 
but we do not yet include any smaller scale time-dependent 
effects (e.g., waves, instabilities, or turbulence) on the inter- 
planetary evolution of these flux tubes. 

For a frame co rotating wi th constant angular velocity SI, 
iWhangld 19801) and lHuld 19931) described the time-steady MHD 
conservation equations in the ecliptic plane (6 = ir/2). In gen- 
eral, we solve for the six magnetofluid parameters p (mass 
density), v r (radial velocity), (azimuthal velocity), P (gas 
pressure), B r (radial field strength), and (azimuthal field 
strength) as a function of r and </>. Mass conservation is given 
by 

M^) + ;|^) = - (7) 

The two nontrivial components of momentum conservation 
are 

dv r Vft, dv r 1 dP GMq 

Vr lT + JH + -7T + ^^ = 
or r o<p p or r z 

B 

-r-MV xB) e + ^ + fi(2vrf, + f7r) (8) 

47T/9 r 

dv^ <9v 1 dP B r v r v^ 
or r Otp pr oq> 4irp r 
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where G is the Newtonian gravitation constant, M Q is the so- 
lar mass, and the latitudinal component of the curl of B is 
given by 

1 dB r OBa, Bj) 
(VxB) e = -— (10) 
r oq> or r 



The internal energy equation for pressure is 

1 <9v0 



dP va, dP 7 
Or r o<p 



— — 

r 2 dr 



= 0. 



(11) 



and the equation of state jP = pa 2 defines the single-fluid 
sound speed a. We choose a constant value for the adiabatic 
exponent 7. The divergence-free constraint on the magnetic 
field gives rise to an effective conservation equation for B r , 
which is 

^(^,) + ^=0. (12) 
r L or r o<p 

Lastly, the "frozen-in" MHD assumption that the magnetic 
field vector remains parallel to the corotating velocity vector 
allows us to write B^, = v^B r jv r . Making use of the defining 
equations for a and B^ leaves 5 equations for 5 unknowns (p, 

V,., V0, P, B r ). 

The MHD equations as described above involve a number 
of key approximations. Our neglect of explicit time depen- 
dence results in the elimination of MHD waves from the mod- 
eled system. However, these equations do contain the terms 
needed to model information pro pagati on on MHD charac- 
teristics, which was used by iHul y993) to successfully pre- 
dict the formation of shocks in the outer heliosphere. Note 
that we also neglect the source terms in the momentum and 
internal energy equations that were required to produce tur- 
bulent coronal heating and wave-pressure acceleration in the 
ZEPHYR models. Thus, we are essentially replacing the ef- 
fects of turbulent heating at heights above our inner boundary 
radius (i.e., r > 2QR Q ) with a simple gas pressure gradient 
determined by the adopted value of the polytropic exponent 
7- 

To solve Equations (ITli-(fT2l. we first constructed the 5x5 
matrices that contain the coeffi cient s multiplying the r and 
<p partial derivatives. Following IHul (11993b . we inverted the 
first of those matrices and solved for each of the r derivatives 
by themselves. To step upward from a lower boundary radius, 
we impleme nted the first - order upwi nd differencing techniqu e 
described bv lPress et all (119921) andlRilev & Lionellol(l2011l) . 
The approach of iRilev & Lionellol (1201 lb was to solve an 
inviscid Burgers' equation that corresponds to Equation dg} 
in the limit of a — > and Va — > 0. We extended their up- 
wind differencing technique to the full set of corotating MHD 
equations, and we used the standard Courant-Friedrichs-Lewy 
(CFL) criterion (treating r as the time-like variable and <fi as 
the space-like variable) to set the radial step size. 

To test our numerical implementation of the MHD equa- 
tions, we solved the example problem of a single "trape- 
zoidal" shaped high-speed stream that was posed by IHul 
(1993). Starting at a lower boundary radius of 0.3 AU, the 
wind speed was doubled (from 300 to 600 km s -1 ) in a 60° 
wide band of longitude, the density was correspondingly re- 
duced such that p oc v~ 1,5 , and the gas pressure was assumed 
to remain constant. This model used an ideal adiabatic value 
of 7 = 5/3. Figure [TTIshows the results of evolving these con- 
ditions from 0.3 to 1 AU, and i t compare s favorably to the 
properties shown in Figure 8 of |Hu| (fl993b . Instead of plot- 
ting the azimuthal velocity in the corotating frame, Figure 
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Figure 11. Azimuthal profiles of plasma properties computed for the IHul 
1 1993) test problem and shown at r = 0.3 AU (blue dotted curves) and r — 1 
AU (black solid curves). Shown are: (a) radial wind speed v r , (b) azimuthal 
speed in the inertial frame u^, (c) proton number density n = p/m p , and (c) 
one-fluid temperature T. The azimuthal angle <j> has been shifted by a con- 
stant offset to have the scale line up with that of Figure 8 of Hu ( 1953). 

[TTTb) shows the azimuthal velocity in the inertial frame, 

W0 = + fir , (13) 

which exhibits departures from zero only because of the 
stream-stream interactions that grow in magnitude with in- 
creasing heliocentric distance. 

Next, we made use of the ZEPHYR solar wind models as 
inner boundary conditions for a two-dimensional model of the 
heliosphere in the ecliptic plane. We focused only on the stan- 
dard parameter choice of T2 = 0.5 MK that was used for the 
models shown in Figure [10] We extracted the density, wind 
speed, gas pressure, and magnetic field strength as a func- 
tion of <j> at an inner boundary radius chosen initially to be 
r = 2QRq. We chose to set = —Or at this radius, which 
is equivalent to assuming no rotation in the inertial frame 
(m^ = 0) as is appropriate for a height well above the solar 
wind's Alfven point (see, e.g.. IWeber & Davisl Il967). Al- 
though we assume time-steady corotation of MHD "patterns" 
in this model system, that is not equivalent to the rigid rotation 
of the fluid itself (i.e., assuming = would be unrealistic). 

The fully resolved SOLIS flux tubes subtend only 80° out 
of a full 360°, so we interpolated the plasma parameters lin- 
early throughout the unused 280° of longitude. We used a 
fine azimuthal grid with 39279 zones and constant separation 
A(p = 0.00917°. We set 7 = 1 .2 to account for ongoing gradual 
heating in the range of distances to be modeled in the eclip- 
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tic plane. This exponent is roughly consistent with an in situ 
temperature dependence of roughly T oc r~ 1//2 , which is close 
to what is seen for the average of pr oton and electron tem per- 
atures in the inner heliosphere (e.g., Cranm er et al.H2009b . 

When running the two-dimensional MHD code for the 
high-resolution set of SOLIS flux tubes, we found that in 
some cases the standard CFL criterion did not govern the nu- 
merical stability of the upwind differencing scheme. In these 
cases, even the use of a radial step size several orders of mag- 
nitude smaller than dictated by CFL would not stabilize the 
system's evolution. We found that Equation (0, the azimuthal 
momentum conservat ion e quatio n, was responsible for this 
instability. Although [Hu| d 19931) found that standard MHD 
characteristics are valid for these equations in the supersonic 
and super-Alfvenic limits, it is possible that the Coriolis and 
curvature terms in Equations d8}-((9} modify the system's ef- 
fective "wave speeds" such that the smaller azimuthal flows 
require special treatment. In practice, we found two alternate 
ways of producing stable models: 

1 . One solution was to assume that deviations from strict 
"Parker spiral" corotation (i.e., = -fir) can be ne- 
glected. Figure [TTT b) showed that departures from 

= at large-scale stream interfaces are typically 
subsonic in magnitude and thus probably unimportant 
to the overall dynamics of CIRs. By making the as- 
sumption that V0 has a known analytic value, we were 
able to exclude Equation (O from the system of equa- 
tions solved by the code. The strict upwind differ enc- 
ing method suggested by Rilev & Lionello (2011) re- 
mained stable for CFL-compliant radial step sizes. We 
refer to this as our "low-diffusion" model. 

2. Another solution was to solve the equation using 
a method that was both more numerically stable and 
more diffusive. W e applied the Lax method (see, e.g., 
iPress et al.l 19921) to the solution of Equation (0, but we 
kept the less diffusive upwind differencing method for 
the other four MHD equations. The inclusion of depar- 
tures from strict corotation (i.e., allowing the system to 
develop longitudinal variations in v^) also gives rise to 
additional physical diffusion in <f>. Both types of diffu- 
sion contribute to the smearing out of sharp structures 
in longitude. Thus, we refer to this case as our "high- 
diffusion" model. 

Below we show results from both models. We believe the 
true solutions to Equations (ITli- (fT2l should lie in between the 
results of the low and high diffusion cases. 

Figure Q~2] shows how the longitudinal profile of radial ve- 
locity v r evolves with increasing radial distance. The "left- 
ward" evolution of structures in (j> is due t o the curvature of 
the ov erall Parker spiral. As predicted by Wan g & Sheeleyl 
(2003) and others, the strongest and narrowest velocity peaks 
at the inner boundary are rapidly smeared out by interactions 
with surrounding slow wind. At the intermediate distance of 
40 Rq shown in Figure [T^ b). the high-diffusion model looks 
somewhat like a smoothed version of the low-diffusion model. 
However, by the time the streams reach 1 AU (Figure Q~2jc)), 
there has been enough dynamical evolution to make the two 
models appear quite different from one another. 

Our use of the polytropic exponent 7=1.2 resulted in some 
extra acceleration of the wind from 20 to 215 Rq. The mean 
wind speed at 1 AU is roughly 10% higher than the speeds 
at the inner boundary. Note that the set of input ZEPHYR 
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Figure 12. Radial wind speeds for MHD models computed from ZEPHYR 
solar wind properties, (a) The inner boundary radius r = 20Rq contains the 
same data for both models. At distances of (b) r = 40i?o and (c) r = 215 Rq = 
1 AU, the high-diffusion model (black curves) and low-diffusion model (red 
curves) differ from one another, (d) In situ data for the modeled time period is 
shown from the multi-source OMNIWeb interface. (One data point is shown 
every 5 minutes, which is downsampled from the original time resolution of 
1 minute.) 



models exhibited a slightly larger degree of acceleration in 
the inner heliosphere, with their original speeds at 1 AU be- 
ing roughly 25% higher than their speeds at 20 Rq. Some of 
the relative lack of acceleration in the two-dimensional MHD 
model could be due to the accumulative loss of kinetic energy 
at multiple colliding CIRs. 

Figure fTZl d) shows the OMNI solar wind speed in the eclip- 
tic plane (measured mostly by ACE/SWEPAM during this 
time period) mapped to longitude using the Carrington rota- 
tion parameters discussed in Section [2] There is not much 
detailed agreement with the MHD models, but the persistent 
local minimum in wind speed centered on <j> s» -40° may be 
a relevant similarity between the data and models. We should 
clarify that the MHD models do not include the full 360° of 
longitude, so there may be large-scale CIR interactions with 
regions outside the high-resolution set of flux tubes that af- 
fect the solutions inside. Also, the agreement between global 
models and in-ecliptic data often depends c rucially on wha t 
was assumed for the Sun's polar fields (e.g., Jian et al. 201 1). 
The magnetograms that are used typically for synoptic recon- 
structions of the global field (e.g., MDI on SOHO, Kitt Peak, 
WSO) often undergo detailed deprojections and high-latitude 
corrections that have not been attempted for the SOLIS data 
used here. 
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Figure 13. Simulated time variability of the low-diffusion (red curves) and 
high-diffusion (black curves) models that would be seen by a spacecraft at 
1 AU. Shown are: (a) radial wind speed v r , (b) azimuthal speed in the iner- 
tial frame Uj,, (c) proton number density n, (d) radial field strength B r , (e) 
azimuthal field strength B^, and (f) Parker spiral angle 8,b- 

Figure [13] shows a summary of other plasma parameters at 
1 AU, here plotted as a function of time as would be mea- 
sured by a stationary spacecraft. We used Equation (|3) to 
translate longitude into time, and for simplicity we normal- 
ized t to zero at the left edge of the plot. Table Q] gives the 
mean values of several relevant parameters at 1 AU for the 
two models. One can see more clearly that the low-diffusion 
model contains sharper features, often with larger variability, 
that tend to be more smeared out in the high-diffusion model. 
The tight correlation between B r and is a preliminary in- 
dication that the dominant type of magnetic variability is in 
the magnitude of the field and not in its direction. Also, the 
modeled fluctuations in gas pressure and magnetic pressure 
are largely in phase with one another, and are not 180° out 
of phase like one would expect in pressure-balance structures 
( Burla ga et aljl990t iTu & Marscri!994t IVasquez & Hollwed 
11999ft . These correlations are discussed further in Section |6) 
The Parker spiral angle is defined as 9,-g = tan" 1 (B^,/B r ). Fig- 
ure lOf f) shows that the pattern of longitudinal variability seen 
for 6 r B does not depend much on whether departures from 
strict corotation are included in (as in the high diffusion 
model) or not (as in the low diffusion model). 

We found it instructive to simulate global time-distance 
images from the MHD models and compare them to simi- 



Table 1 

Modeled MHD Parameters at 1 AU 



Parameter 


High- diffusion Model 


Low-diffusion Model 


Mean v r (km s ) 


474.3 


481.2 


Mean (km s _1 ) 


^00.2 


-399.0 


Mean n (cm -3 ) 


3.854 


3.656 


Mean |B| (nT) 


2.169 


2.072 


Mean /3 (iWiW) 


7.383 


8.223 


Ek.//E„ 


1.402 


1.612 


Ek.4>/E„ 


0.1407 





Eb,i/E„ 


0.06527 


0.09566 


Eb.^/E,, 


0.04293 


0.06724 


E t h/E„ 


0.5919 


1.131 


Total E/E„ 


2.243 


2.906 
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Figure 14. Simulated time-distance variations of density fluctuations from 
the high-diffusion MHD model. Th e no rmalization of the horizontal time 
axis is the same as that used in Figure [131 Grayscale levels are defined inde- 
pendently at each radial distance to be proportional to relative electron density 
variations. 

lar images constructed from observations of white-light he- 
liospheric imagers. For example, the Sun Earth Connec- 
tion Coronal and Heliospheric Investigation (SECCHI) on 
STEREO contains imagers that probe el ectron density flu ctu- 
ations at distances of order 15-300 Rq (Evle s et alJ2009l) . By 
stacking up radially oriented image slices taken over a range 
of times, one can visualize the mutual interactions of fast and 
slow streams, coronal mas s ejections (CMEs), and shocks in a 
clear way (see, e.g.. | Davies et al. 2009; Rouillard et al.H20lot 
Howard & DeForest 2012). Figure [T4lshows the result, using 
the same time coordinate as in Figure [T3] and using the mass 
density p as the plotted grayscale quantity. At each height, 
the white and black limits of the grayscale map were redefined 
based on the local minimum and maximum density in <f>. Both 
the modeled and observed images show that fine-scale flux 
tube variations in the innermost heliosphere become smeared 
out by CIR-like stream interactions at larger distances. 

Many of the results discussed above are independent of the 
adopted value of the inner boundary distance. As discussed 
above, we chose 20 Rq as a standard baseline. We wanted a 
value low enough that substantial stream interactions would 
not have already occurred, but high enough that most of the 
solar wind's acceleration as modeled by ZEPHYR would be 
included. We also wanted to avoid the sub-Alfvenic region 
in which is undergoing a transi tion from rigid rotatio n to 
angular momentum conservation (Weber & Davis 1967J). In 
order to explore the sensitivity of the parameters at 1 AU to 
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the boundary distance, we ran two other low-diffusion models 
with inner radii of 15 Rq and 30 Rq. Because we did not 
fine-tune the polytropic 7 to reproduce the ZEPHYR code's 
wind acceleration profile at large distances, we found that the 
mean wind speed v r at 1 AU does depend on the choice of 
inner radius. However, in all three models, the mean mass 
flux pv r and the mean radial field strength B r at 1 AU were 
nearly invariant (i.e., they differed by less than 1% from one 
another). Also, the relative fluctuations at 1 AU, measured by 
the ratios of standard deviations to mean values (for both p 
and v r ) were nearly invariant across the three models. 

6. STATISTICS OF MODELED MHD FLUCTUATIONS 

Because the two-dimensional models described above do 
not include any explicit time variability, the resulting "fluctu- 
ations" (as seen in, e.g., Figure [T3l should not be expected 
to resemble propagating MHD waves or turbulent eddies. 
Nonetheless, it is useful to first examine them in the context 
of linear MHD wave theory. We determined the mean energy 
density components of the fluctuations at 1 AU, 



Po(vf) 



(Bf) 



2po 



(14) 



where quantities with subscript refer to mean values taken 
over the 80° of fully resolved longitudes, quantities given as 
(f 2 ) refer to the variance of / over the same range of lon- 
gitudes, and the subscript i refers to either the r or <fi vector 
components. The three terms above refer to kinetic, mag- 
netic, and thermal fluctuations, respectively. Alfvenic fluc- 
tuations would be characterized by equipartition between the 
kinetic and magnetic components perpendicular to the back- 
ground field. Magnetosonic waves tend to have half of their 
fluctuation energy in kinetic form and the other half divided 
between the magnetic and therma l terms (see Figure 2 of 
ICranmer & van Ballegooiien 2012). Pressure-balance struc- 
tures (PBSs) are largely static features that advect with the 
solar wind and carry along variations in gas pressure and mag- 
netic pressure that are in rough equipartition between the Er 
and E± terms given above. 

Table Q] gives the values of the fluctuation energy density 
components for the two models discussed in Section [5] Each 
energy density has been normalized by dividing by a represen- 
tative 1 AU background magnetic energy density E„ = B 2 /%n, 
where B„ = 2x 10~ 5 G = 2 nT. This value of B„ was chosen to 
be of the same order of magnitude as the mean magnetic field 
strength at 1 AU in the models. TableQ]shows that both mod- 
els had significantly more than half of their total energy in the 
form of kinetic fluctuations, which is inconsistent with both 
linear MHD waves and PBSs. Still, the models had an ap- 
proximately magnetosonic balance between the relative rms 
variations in magnetic field amplitude {{B 2 ) x l 2 /Bq) and den- 
sity ((p 2 ) 1 / 2 1 po), with both ratios being roughly equal to 0.30 
for the high-diffusion model and to 0.39 for the low-diffusion 
model. This may be a natural consequence of the frozen-in 
MHD assumption, since the spacing of field lines undergoes 
similar compressions and rarefactions as the gas. 

The total modeled energy density in fluctuations was be- 
tween 2 and 3 times the background value E n . This can be 
compared with measured plasma fluctuations at 1 AU. The 
total energy density of reported MHD turbulence is roughly 
about 3-5 E„ (e.gjTu & Marschlfl99lt i Goldstein et al| [l99l 
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iBruno & Carbonel2005l) . depending on the type of solar wind 
stream. This is in approximate agreement with the models, but 



Figure 15. Radial dependence of magnetic field variance in the inner helio- 
sphere. For the low-diffusion model (red solid curve) and the high-diffusion 
model (black solid curve), plotted values of the component variance 
are nearly identical to the magnitude variance a 2 b (not shown). These are 
compared with Helios measurements of cr? (dashed curves) and ai (dotted 
curves) for data averaged over 12 hour bins (blue upper curves) and 3 hour 
bins (green lower curves). 

we will see below that there are key differences between the 
detailed power spectra of the models and the observations. We 
also found that larger-scale fluctuations in the ecliptic plane 
contribute to a much higher energy density than do just the 
small-scale turbulent motions. For example, over the 16 days 
of OMNI data collected for Figure [T2f d). the total energy 
density computed from all components of Equation (TBI was 
found to be 42.4 E n . Much of this (28.3 E„) was in the kinetic 
energy variance due to alternating fast and slow streams. It is 
clear that the large-scale models presented here do not come 
close to reproducing the full range of MHD variability in the 
solar wind. 

We can leam more about how the models may reproduce 
some aspects of the actual solar wind by examining the ra- 
dial evolution of the magnetic field fluctuations. Figure Q3] 
shows the radial dependence of the "component variance" a 2 , 
which is the sum of the variances of the individual vector com- 
ponents of B. For both the low-diffusion and high-diffusion 
models, this quantity is nearly identical to (i.e., within 5% of) 
the magnitude variance <j\, which is the variance of the scalar 
field strength |B|. We do not plot a\ for the models since the 
curves would be indistinguishable from the a 2 curves. 

Figure Q3] also compares the model results with measure- 
ments made by Helios in the slow wind (v < 500 km s" 1 ). We 
show measured ranges for a 2 and a\ given by Maria ni et al.l 
(1978) for two different methods of binning the data. In the 
real solar wind, a 2 c la 2 b ~ 10, which indicates that the direc- 
tional variations in the vector magnetic field greatly exceed 
the magnitude variations. Although our model does not con- 
tain the small-scale MHD waves that contribute to the ob- 
served a 2 , it does seem to reproduce the large-scale mag- 
nitude variations that contribute to a 2 b . Thus, we make the 
conjecture that the low-frequency variability of magnetic field 
magnitude in the solar wind may be caused by the longitudi- 
nal interaction and evolution of coronal flux tubes. 
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Figure 16. Fourier power spectra of v r fluctuations in the time-steady MHD 
models. The raw FFT result for the inner boundary condition at 20 i?0 is 
shown (light gray curve) along with its 9-point boxcar smoothing average 
(blue curve). Corresponding smoothed spectra at 1 AU are shown for the 
low-diffusion model (red curve) and high-diffusion model (black curve). 

Additional information about the nature of the modeled 
fluctuations can be extracted from their power spectra. Con- 
verting the longitudinal coordinate <j> into time using Equation 
(01, we took fast Fourier transforms (FFTs) of various MHD 
quantities over the 80° of fully resolved longitude and plot- 
ted the power as a function of spacecraft-frame frequency. 
For example, Figure [16] shows power spectra of the fluctua- 
tions of radial velocity v r at both the inner boundary radius 
of 20 Rq and at 1 AU. The curve at 20 Rq is the same for 
the low-diffusion and high-diffusion models because the same 
ZEPHYR boundary condition was used for both. Despite the 
high resolution of the models, the numerically determined 
power spectra were extremely noisy. Thus, for clarity we plot 
9-point boxcar-smoothed curves for all modeled spectra, and 
we show the full, unsmoothed spectrum for only one case: the 
inner boundary spectrum at 20 Rq. The smoothing does not 
alter the overall shape or magnitude of the spectrum. 

The velocity power spectrum at 20 Rq is roughly propor- 
tional to f~ l for frequencies between about 5 x 10~ 6 and 10~ 4 
Hz. This is reminiscent of — but probably not causally related 
to — the spectrum seen for some quantities at frequencies 
below IP" 4 Hz in the solar wind a t 1 AU (Bavassan o et alj 
[1981 IMatthaeus & Goldstein! 1 198ft iNicol et alJI2009l) . The 
modeled spectrum steepens to roughly f~ 2 at higher frequen- 
cies, and the steepening continues as the spectrum evolves 
outward in heliocentric distance. The high-diffusion model 
steepens more than the low-diffusion model. This is expected 
since the former model undergoes more longitudinal smear- 
ing of structure at small spatial scales (i.e., at high spacecraft- 
frame frequencies). The power at the very lowest frequencies 
(f < 3 x 10~ 6 Hz, i.e., timescales greater than a few days) ap- 
pears to grow with increasing heliocentric distance. Although 
there is no formal "inverse cascade" in this model, the trans- 
fer of energy from small to large scales can be the result of 
CIR-like stream collisions that eventually produce merged in- 
teraction regions (MIRs; e.g., Burl aga et al.l l2003). 
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Figure 17. Fourier power spectra of magnetic field fluctuations in the high- 
diffusion (a) and low-diffusion (b) models at 1 AU. Simulated power spectra 
(solid curves) are compared with measured ranges of vector component fluc- 
tuations (blue dotted curves) and scalar magnitude fluctuations (green dotted 
curves). 

Figure [17] shows the modeled magnetic fluctuation power 
spectra at 1 AU without any boxcar smoothing. For clarity, 
the results of the low-diffusion and high-diffusion models are 
shown in separate panels. The gray regions give the ranges of 
measured magnetic power spectra for the slow solar wind in 
the ecliptic plane. The spectrum labeled P c (f) gives the "full" 
MHD turbulence spectrum that is analogous to the component 
variance aj discuss ed above (see, e.g . . iTu & Marschlll995l: 
MatsuietaJJ l2002t iBruno & Carbond 120051: IVasquez et al I 
2007). The spectrum labeled P;,(f ) is the spectrum of fluc- 
tuations in the scalar field strength (Burlaga et al. 1989]), and 
it is analogous to the magnitude variance a\. The agreement 
between the low-diffusion model spectrum and the measured 
range of Pbif) is quite good. This appears to reinforce our 
earlier suggestion that the proposed model of longitudinally 
evolving coronal flux tubes can explain the measured low- 
frequency fluctuations of magnetic field magnitude in the so- 
lar wind. 

It is also evident from Figure [T7]that the low-frequency 
part of the measured fluctuation spectrum P c (f) contains sig- 
nificantly more power than is found in the modeled spectra at 
1 AU. The f~ l fluctuations are o ften claimed t o be fossil rem- 
nants of coronal variations (e.g., Roberts 2010). However, we 
can tentatively conclude from these simulations that if this is 
true, it is more likely for them to represent propagating waves 
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and not just the "passive" survival of coronal flux tubes. 



7. DISCUSSION AND CONCLUSIONS 

The aim of this paper was to begin exploring how MHD 
fluctuations in interplanetary space are related to the well- 
known filamentary flux-tube structure of the solar corona. 
Starting with a time-steady potential-field model of the QS 
corona allowed us to model the properties of a "bundle" of 
open field lines with footpoint separations that extend down 
to scales smaller than those of both granules and intergran- 
ular lanes. This let us study magnetic field variations on 
spatial scales that are still unresolvable with state-of-the-art 
MHD simulations of the global corona and heliosphere. After 
adjusting the field strengths found from the PFSS approxi- 
mation, we produced one-dimensional models of the coronal 
heating and solar wind acceleration along these field lines. 
These models assumed that the energization of the corona 
comes from the dissipation of MHD turbulence, but in prin- 
ciple any model of plasma heating could be applied. We also 
used the one-dimensional results as lower boundary condi- 
tions to a model of radial and longitudinal heliospheric evolu- 
tion, and we predicted the time-steady formation of CIR struc- 
ture from the coronal flux tubes. The resulting fluctuations 
in magnetic field magnitude greatly resemble corresponding 
fluctuations that have been measured in the low-latitude slow 
solar wind. 

One original goal of this paper was to study the ori- 
gin of low-frequency MHD fluctuations in the solar wind. 
Power spectra obeying a relatively flat f frequency depen- 
dence are often suggested to be the turbulent "energy con- 
taining range" and of direct solar origin, but their cause is 
still be ing debated (M atthaeus & Goldstein 1986 : IVelli etalJ 
| 1989t iMatthaeus et all 120071: iDmitruk & Matthaeus! 120071: 
IPmitruk et al.ll201 It IVerdini et al.ll2012bl) . There is also still 
disagreement over whether the statistical propertie s of MHD 
disco ntinuities are due to flux tubes of solar origin (iBorovskvl 
or are the by-product of an ongoing turbulent cascade 
dGreco et al.ll2.008n . Our ecliptic-plane models were not set 
up to reproduce inherently time-variable fluctuations, but they 
do help put limits on the ability of small-scale coronal varia- 
tions to survive the effects of stream-stream CIR-type interac- 
tions. They also act as a basal level of expected background 
fluctuations to which must be added any propagating waves 
or turbulent motions. Interactions between waves and static 
structures have been invoked as pos sible means of mode con- 
version and coron al heating (e.g., iGogoberidze et all 120071 
IE vans et all 1201 2|) . so knowing more about the background 
variations is important to help constrain those theories. 

As emphasized in Section[2l our time-steady description of 
the coronal magnetic field was meant only to be a snapshot 
that is representative of the hi gh-resolution topology th at ex- 
ists at any one time (see also Jendersie & Peter! 120061) . The 
validity of this kind of model depends crucially on the scales 
being resolved. Because the smallest structures will have the 
shortest lifetimes, it becomes increasingly necessary, when 
going to smaller scales, to ta ke into account the evolving 
nature of the magnetic carpet. Cranm er & van B allegooiien 



3 Above a certain height, the magnetic carpet's evolution may. be slqw 
enough that a parcel of solar wind can accelerate up to that height without its 

footpoints being appreciably deformed. In that case, our "snapshot" method 



(2010), iMever et all d2012l) . and others presented Monte 
Carlo simulations that may be useful to "connect" to the so- 
lar wind acceleration models of this paper. At the very least, 
such simulations can be used to provide insight about how 
rapidly the field geometry evolves as a function of heightQ 
Of course, fully time-dependent and multidimensional tech- 
niques will be needed to describe the self-consistent interac- 
tions between flux tubes, waves, turbulent fluctuations, and 
fast/slow streams in the heliosphere. 

We also note that our assumption of a potential field for ex- 
trapolating the measured photospheric magnetogram data into 
the corona is likely to be inaccurate in many key ways. Us- 
ing an improved MHD treatment could naturally incorporate 
many of the non-potential effects that we attempted to account 
for in our Bi enhancement term. In addition, the PFSS tech- 
nique is known to eliminate nearly all of the smallest- scale 
longitudinal flux-tube structure in the field at and above the 
source surface. In our case, this resulted in an ina bility to test 
the alt erna te flux- t ube e xpansion formalism of Fuii ki et all 
d2005l) and lSuzukil (120061) . who suggested that B ss oc B base /f ss 
may be a better quantity to correlate with the wind speed than 
/ ss itself. 

Despite the above warnings, we believe the high-resolution 
magnetic field models constructed for this paper can be useful 
tools for studying the statistical properties of spatial plasma 
fluctuations that exist sufficiently near the Sun that other ef- 
fects (e.g., stream interactions or turbulence) do not smear 
them out. These variations can be compared straightforwardly 
to existing remote-sensing observations. For example, we 
would like to know if these inter-flux-tube variations con- 
tribute to the m easured nonth ermal widths of coronal emis- 
sion lines (e.g.. iGuptal 120 1 dl) or to the lar ge velocity vari- 
ances obtained from radio scintillation data (Har mon & Coles! 
2005). Both of these are commonly interpreted as being 
caused by the line-of-sight integration through a corona filled 
by Alfven waves, but the effects of "static" flux tubes re- 
main to be disentangled. We would also like to simulate 
the MHD fluctuations at the intermediate distances to be 
probed by Solar Probe Plus (r > 0.05 AU) and Solar Orbiter 
(r > 0.28 AU). This requires a more self-consistent calcu- 
lation of to model the turnover fr om rigid corotation to 
angular momentum cons ervation (e.g., Web er & Davisl[l967[ 
iPriest & Pneumanl[l9 74), and it may also require including 
"true" heating and acceleration source terms into the two- 
dimensional code described in Section|5] 
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APPENDIX 

POTENTIAL FIELD SOURCE SURFACE MODEL 

In this paper we compute the magnetic field assuming the corona is current-free (i.e., a potential field). We use a spherical 
domain extending from the solar surface up to a radial distance R ss , where the field is assumed to be ra dial. Low resolution 
versions of such potential field source surface (PFSS ) models have been used extensively in the past (e.g..lSchatten et al.l 119691: 
Altschuler & Newkirk|[l969l iWang & Sheelevlll990t lArge & Pizzdl2000t iLuhmann et allllOOl IWang & Sheelevll2006l) . In the 
present work we need high spatial resolution in the low corona to accurately trace the roots of the open field lines on scales 
below ~1 Mm in the photosphere. Using the standard method, which is based on spherical harmonic expansion, very large grids 
covering the entire Sun would be required. However, high resolution is needed only in a limited region on the solar surface. 
Therefore, we use a different computational method i n which the high resolution grid extends only over some limited area. 
The basic methodology is described in Appendix B of van Ballegooijen et al. (2000), and can be summarized as follows. The 
magnetic field B(r) is written in terms of a vector potential A(r), which is given by 

1 <9<J/ 

A r = 0, A g = ; — — — , (Al) 
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(A2) 



where (r,9,(j)) are spherical coordinates, Bo is the monopole component of the imposed radial field, and ^>(r,9,(f>) is a scalar 
function that must satisfy the following partial differential equation: 

+ Vi* = 0. (A3) 

The gradient operator Vj_ operates only on the coordinates perpendicular to the radial direction. Thus, the potential field is given 
by 

B(r) = £o(^)Vv(f). (A4) 

At the longitudinal boundaries of the high-resolution domain, B$ either vanishes or is p eriod ic in longitude; at the latitudinal 
boundaries we require Bg = 0. Due to these boundary conditions, the solution of Equation (IA3t can be written as a superposition 
of discrete eigenmodes, but the eigenmodes are not spherical harmonics and must be computed numerically. The mode amplitudes 
are chosen such that B(r) matches the imposed radial magnetic field at the solar surface (r = Rq) and satisfies d 1 ^ jdr = at the 
source surface (r = R ss ). 

Two important improvements have been made to the basic method. First, we use a variable grid. At certain heights in the 
corona the grid spacing is doubled in all three directions, all owing us to cover t he full height range of the model with a relatively 
small number of grid points (for details see Appendix A of Bo bra et aLl l2008). Second, a description of t he global coro na has 
been added to improve the side boundary conditions on the high-resolution domain (see Section 4.1 of ISu et ail 1201 ll) . The 
high-resolution and global models are computed in such a way that the magnetic field is continuous across the side boundaries 
of the high-resolution domain. Hence, field lines can be traced through the side boundaries without any unphysical kinks in the 
field lines. 

As described above, we used a SOLIS line-of-sight (LOS) photospheric magnetogram taken on 2003 September 4 at 16:30 
UT, as well as the SOLIS synoptic map for the relevant Carrington Rotation (CR 2007). The LOS magnetogram was converted 
into a longitude-latitude map of the radial field B r (R Q ,9, <fi) in the high resolution part of the grid. This conversion assumes that 
the magnetic field on the Sun is nearly radial, which is a good approximation for the Quiet Sun. The high resolution grid has 
1024 x 512 cells on the photosphere, and the grid spacing is 0.0018cos XRq, where A is the latitude. The source surface is located 
at R^ = 2.46 Rq. The synoptic map from CR 2007 was used for constructing the global part of the PFSS model. 
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